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We study the linear instability with respect to necking of a filament of polymeric fluid undergo- 
ing uniaxial extension. Contrary to the widely discussed Considere criterion, we flnd the onset of 
instability to relate closely to the onset of downward curvature in the time (and so strain) evolution 
of the zz component of the molecular strain, for extension along the z axis. In establishing this 
result numerically across five of the most widely used models of polymer rheology, and by analytical 
calculation, we argue it to apply generically. Particularly emphasized is the importance of polymer 
chain stretching in partially mitigating necking. We comment finally on the relationship between 
necking and the shape of the underlying steady state constitutive curve for homogeneous extension. 



Understanding the rheology (flow properties) of poly- 
meric fluids is central to their processing and perfor- 
mance. Many commercially important flows are dom- 
inated by extensional components: fibre spinning, film 
blowing, extrusion, and ink jet printing provide good 
examples. In fluid dynamical terms, extensional flows 
cause material elements to separate exponentially quickly 
and so subject the underlying macromolecules to extreme 
stretching and reorientation. They are thus highly sensi- 
tive to underlying molecular details: linear vs branched 
polymer chains, for example. Indeed many nonlinear flow 
features arise only in extension, which thus provides a 
crucial benchmark for theories of polymer rheology. 

In many polymeric fluids, a state of uniform flow be- 
comes unstable when the flow rate exceeds the rate 1/t 



on which the underlying molecular structure relaxes 
In flows dominated by shear, for example, the phe- 
nomenon of shear banding Q arises widely for shear rates 
7T > 0(1). In a linear stability analysis, the criterion for 
onset of a shear banding instability is (usually) that the 
shear stress is a decreasing function of shear rate. 

In extensional rheology [1,|3| a common protocol is that 
of filament stretching [5| , in which a cylindrical sample of 
fiuid is drawn out in length. Though the aim is to achieve 
uniform extensional fiow for benchmarking against the- 
ory, complications often arise. Particularly serious is the 
widespread observation of necking instabilities (Fig. [T|), 
in both experiment and simulation [lj-ll7|. These 

lead to heterogeneous deformation and even failure at 
only modest strains: any small indentations in cross- 
sectional area become ever more pronounced until the 
sample breaks. This hinders attempts to characterise 
these fiuids scientifically, and process them commercially. 

Although a long standing problem 



18H291, necking re- 



mains poorly understood. Crucially lacking is any reli- 
able criterion for its onset, of the same stature as that 
given above for shear banding. Popularly discussed is the 
Considere criterion [s^l , which predicts necking if the ten- 
sile force decreases with extensional strain: dF/de < 0. 
But taking F to depend only on e in this way assumes the 
flow purely elastic: it cannot account for the dependence 
of necking on the strain rate e in these viscoelastic fluids; 



nor predict the rate at which necking sets in. 

Here we show that the Considere criterion in fact does 
not apply in most regimes of polymeric flow. By a linear 
stability analysis, we demonstrate the onset of necking in- 
stead to relate closely to that of downward curvature in 
the time (or equivalently strain) evolution of the zz com- 
ponent of the molecular conformation tensor, for stretch- 
ing at constant rate e along the z axis. In most regimes, 
this further corresponds to downward curvature in the 
evolution of the tensile stress. In establishing this result 
numerically across flve of the most widely used models of 
polymer rheology, and by analytical calculation, we argue 
it to apply generically. The other central contribution of 
this Letter is to demonstrate the crucial role played by 
polymeric chain stretching in partly mitigating necking. 

We consider a polymeric fluid in inertialess flow, obey- 
ing the force balance condition = V.T = V.{'S + 2riD — 
PI). The total stress T has a viscoelastic contribution 
S from the polymer; a Newtonian contribution of vis- 
cosity r] from the solvent; and an isotropic pressure held 
P set by mass balance V.v — for incompressible flow. 
Here v is the velocity held and D the symmetric part 
of the velocity gradient tensor = daVp. We assume 
the polymeric stress S = GQ[k{W)W — I] with Go a 
modulus and W a dimensionless tensor characterising 
the conformation of the chainlike polymer molecules. (In 
the simplest cartoon these are taken as dumbbells with 
span R and W oc (RR).) This has generalised dynamics 

{dt+v ■V)W ■ K + k'^ -W - -R{W), (1) 

r 

representing a competition between driving out of equi- 
librium by flow (terms in k) and relaxation back to equi- 
librium W = J on a timescale r. 



FIG. 1: Cartoon of necking. 
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FIG. 2: a) Homogeneous constitutive curves for steady state 
uniaxial extension. Clockwise: Oldroyd B; Giesekus {a — 
0.001); fene-P {S = 0.001); stretch Rolie-Poly (tr = 0.0316) 
with zero extensibility; the same with finite extensibility (/ = 
10~*); non-stretch Rolie-Poly. (Curves for Giesekus and fene- 
P indistinguishable.) b) Tensile stress cte vs time in startup of 
uniaxial extension in the stretch Rolie-Poly model with tr — 
0.0316,/ = 10"'' for i = 10", n = -1,-1/2,0,1/2,1,3/2,2 
(curves upwards). Inset: corresponding tensile force. 



We perform calculations for five concrete functional 
choices for k{W) and R{W), each correspondin g t o 
a particular model used widely in the literature [3lj. 
The Oldroyd B model has fc = 1, = W - J. 
The Giesekus model differs from this in having 
R^{W -I) + a(W - If with < a < 1; the Fene-P 
(in which we assume small 5) in having R — k{W)W — I 
with A; = 1/(1 - 5T) and trace T = W^^ -t- Wyy + W^^; 
the Rolie-Poly ^ R = W ~ I + 2{l ~ 
.JijT) [W + /?(T/3)^-p(W - /)] r/ [rR,(l - /r/3)]. 
Finite tr in the Rolie-Poly model allows stretching of 
the underlying polymer chains by the flow field. Indeed 
for eTK > 1 infinite stretch can develop if / = 0; / > 
restores finite stretch; tr — > disallows it entirely. 

We use units in which t ~ 1, Go — 1. This leaves as 
parameters the Newtonian viscosity rj, which we take to 
zero; and a (Giesekus), 5 (Fene-P), /3,tr , ^rp,/ (Rolie- 
Poly) in which we use /3 = 0.0, Jrp = — 1/2 following [32i| . 

We consider a long slender [33[ tube of fluid subject 
to uniaxial extension along the z axis. At any time 
t it has cross sectional area profile A^z^t) and area- 
averaged fluid velocity V{z,t) in the z direction, with 
extension rate e(z, t) — dzV{z, t). The mass balance con- 
dition is then dtA{z,t) — —dz(AV), and force balance 
F{t) = A[k{T){Z - X) + 37je] = Aa^ (neglecting sur- 
face tension). This defines the tensile stress cteI-z, t) given 
force F{t), which is uniform along the filament. Z — Wzz 



and X = hWr 



(so Z + 2X = T) obey 



dtZ{z,t) + VdzZ = 2eZ - R{Z,T)/t, 
dtX{z,t) + VdzX = -eX - R{X,T)/t, (2) 

with R{Z,T) = Z - 1 (Oldroyd B); i?(Z,T) 
Z - I + a{Z - 1)2 (Giesekus); R{Z,T) 
k{T)Z - 1 (Fene-P); R{Z,T) = Z - 1 + 2{l - 

7|) [z + /3(f)*-- (Z - 1)] /rR(l - ff) (Rolie-Poly). 



Given a time-dependent separation of the sample ends 
that imposes a global strain e(t) averaged along the fil- 
ament, we transform to the affinely co-extending, co- 
thinning frame by defining u = ze~'^^*\ v{u,t) — 
y(z,t)e-'^(*) and a{u,t) = A{z,t)e^'^*\ We then have 
mass balance, force balance, and viscoelastic dynamics: 
dta = -du [{v - eu)a\ , (3) 
= duHkiT){Z -X) + 377e)] = du[aaE], (4) 
dtZ = -{v ~lu)duZ + 2eZ - R{Z,T)/t, 
dtX ^ -{v ~eu)duX ^ eX - R{X,T)/t. (5) 

For a constant global extension rate e commenced at 
time t — 0, a homogeneous "base state" in which uniaxial 
extension is (artificially) maintained uniformly along the 
filament, regardless of whether it is in practice unstable 
to heterogeneous necking, is prescribed by e{u,t) = e, 
v{u,t) = ue, a{u,t) = ag = 1, and the homogeneous 
solutions Z{t),X{t) of Eqns. Elgiven Z(0) = X(0) = 1. 

In any regime of finite extensional viscosity this base 
state attains as t — >■ oo a steady state in the co-thinning, 
co-extending frame. (In the laboratory frame the sample 
exponentially extends and thins.) This is described by 
a homogeneous constitutive curve crE(e) of tensile stress 
vs strain rate. Fig. [5^1- (For this uniform base state we use 
symbols e and e interchangeably.) The Oldroyd B model 
has divergent viscosity for e — ^ 1/2. This divergence is 
avoided (narrowly, for small a, 5) in Giesekus and Fene- 
P (which reduce to Oldroyd B for a = 0, 5 = 0). The 
stretch Rolie-Poly model with / = likewise has diver- 
gent viscosity for etr — 1, avoided narrowly for small 
/ > and entirely in the non-stretch model tr =0. 

The time-dependence of the tensile stress ue in this 
homogeneous base state as it evolves towards the steady 
state just described is shown in Fig. [21d for the stretch 
Rolie-Poly model with finite extensibility. The corre- 
sponding force F{t) — a-E{t)A{t) (inset) displays a max- 
imum due to the rise in (Te after the inception of flow 
being later overcome by exponentially declining area. 

A maximum in force vs time t directly also signifies a 
maximum vs strain e = it, given constant e. For a non- 
linear elastic solid, the Considere criterion decrees such 
a maximum to herald departure from uniform extension 
and onset of necking. We now examine whether this cri- 
terion also applies to polymeric fluids, as often suggested. 

To do so, we analyse the linear stability of the time- 
dependent homogeneous uniaxially extending base state 
just described by linearising in small heterogeneous fluc- 
tuations '^^[5aq{t),5eq{t),5Zq{t),5Xq(t)]eyi]}(iqu) about 
it. To flrst order in the amplitude of these we have 
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FIG. 3: Contour map for growth of fluctuations in area: 
5a{t) / 5a{G) — 2" for n — 1, 2, 3 • • • sliown by thin solid lines 
from left to right. Dotted lines, open circles, thick solid lines 
show location respectively of f = 0, ct'e ~ 0, Z = 0. a) 
Oldroyd B, b) Giesekus {a = 0.001), c) Fene-P (S = 0.001). 
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FIG. 4; Contour maps 5a{t)/5a{0) = 2" for n = 1, 2, 3 • • • , 
thin solid lines left to right. Dotted lines, open circles, thick 
solid lines show location of / = 0, ite = 0, Z = 0. a) Non- 
stretch Rohe-Poly, b) Stretch Rolie-Poly (tr = 0.0316, / = 
0.0), c) Stretch Rolie-Poly with finite extensibility, / — 10~*. 
Red line in a) shows analytical prediction for divergence in 
Sa. 



will grow (or decay) equally, though we return below to 
discuss what shape can be expected in practice. 

If the time-dependence of the base state were disre- 
garded, instability to necking would correspond to any 
eigenvalue of this matrix having positive real part. How- 
ever its entries Z, X,a-E,T refer to the time-evolving base 
state, rendering the eigenvalues time-dependent: indeed, 
these can start out all negative before one later goes pos- 
itive. Numerics and analytics (not shown) reveal this 
sign change to correspond to excellent approximation to 
onset of downward curvature Z < in the time (and so 
strain) evolution of the base state Z{t) and also (apart 
from the small discrepancy seen in Fig. [3]) of the tensile 
stress, (t'e < 0. This strongly suggests that necking will 
arise in any regime oi Z < 0,crE < 0. Some caution 
is needed, however, because the eigenvectors also evolve. 
To examine rigorously the onset of necking, therefore, 
we explicitly integrate the linearised equations numeri- 
cally. (Nonlinear effects must eventually become impor- 
tant, but when depends on the size of the seeding pertur- 
bation, which is not specified.) Resulting contour maps 
of Sa{t) /Sa{0) are shown in Figs. |3l IH with the criteria 
Z = 0,(Te = 0, F = for comparison. 

In the Oldroyd B model. Fig. [3^, 6a{t) grows without 
bound for e < 1/2, albeit slowly, predicting that necking 
should eventually arise in any experiment of long enough 
duration, even causing filament failure unless mitigated 
by nonlinear effects. For e > 1/2, in contrast, Sa(t) grows 
only weakly before quickly saturating, signifying stability 
against necking. This is consistent with the divergence in 
the underlying constitutive curve: for e > 1/2 the base 
state Z{t) and (TE(i) diverge in time with ever upward 
curvature and the eigenvalue remains negative. (For e < 
1/2, Z{t) and (Te(0 curve down towards steady state.) 

The Giesekus and fene-P models avoid this viscosity 
divergence: in the evolution of the base state an early 
time regime of Oldroyd B-like upward curvature (Te > 
0, Z > 0, inside the nose-shaped region in Figs. |3)3,c. 



later gives way to downward curvature toward a high 
viscosity steady state. Surprisingly Sa does grow inside 
the nose, even though Z > and all eigenvalues are neg- 
ative here. This growth is however weak enough to go 
unnoticed, assuming the perturbation that initially seeds 
the instability to be small. Rapid exponential growth 
and observable necking set in only after the nose, where 
Z < 0, (Te < and an eigenvalue is positive. 

Can these contour maps be understood analytically? 
Given an evolution equation dtCTE — iG{aE,T) — 
H{a-E,T)/T for the base state's tensile stress, one can 
show (approximately) for Oldroyd B and Giesekus that 

Sa{t)/Sa{0) ~ f{x) with / (1 + hxf^^ and x — a^loE- 
Here h — H/ tcte — e varies only mildly compared to x in 
Giesekus, and is constant for Oldroyd B. In any regime 
of finite viscosity (Oldroyd B for e < 1/2 and Giesekus 
for all e) the argument x — ^ oo as t — > cx) due to the 
downward curvature to steady state ct'e ^ 0. Combined 
with J{x) — >■ oo as a; — 7- oo, this gives (5a — > oo as t — >■ oo: 
necking fluctuations grow without bound, as seen numer- 
ically. In contrast the divergent viscosity of Oldroyd B 
for e > 1/2 gives x — > const, as t — >■ oo and 5a saturates. 

We now use these findings in the above phenomenolog- 
ical (but widely used) models as pedagogical backdrop to 
the more realistic Rolie-Poly model, focusing particularly 
on the dramatic role of polymeric chain stretch. 

Without chain stretch, tr — 0, the result 5a{t) = 

5a{Q)f{x) with / — {\ A-hxY^^ , x — cte/ct'e and b — 
l/r — e holds exactly. Finite extensional viscosity fur- 
ther gives X — 7- oo as t — > cxD for all e, and so unbounded 
growth of 5a. In fact for e > l/r, / diverges at finite x 
giving a dramatic finite time divergence in 5a. An analyt- 
ical prediction for this is shown in Fig. |4^; the numerical 
contours indeed accumulate at it. Clearly, suppression of 
chain stretch confers spectacular necking instability. 

For finite tr the base state shows divergent extensional 
viscosity for e — > 1/tr,, associated with the development 
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of infinite chain stretch. This is reminiscent of Oldroyd 
B for e — 1/2. Correspondingly, necking is suppressed 
altogether for e > 1/tr (Fig. |3)d). Even for e < I/tr, 
the finite-time divergence of Fig. disappears. Chain 
stretch thus strongly mitigates necking. 

Although instructive, infinite chain stretching is clearly 
unphysical so we now introduce a parameter / to cure 
it [32], analogous to 6 in fene-P. With this, the regime 
of upward curvature in the base state, Z > 0, cte > 0, 
is again confined to a nose-shaped region (Fig. \^). As 
in Giesekus and fene-P, Sa grows only slowly inside this: 
exponential growth sets in only beyond it, once Z < 0, 
(Te < 0. The nose's tip at ctr = 0(1) gives a maximum 
in onset strain vs. e that closely resembles experimental 
data [1, , as collated onto a master plot in Ref . 01 • 

Taken together. Figs. [3] and 2] tell us how heterogeneity 
grows in a necking instability intrinsic to the material's 
rheology. To specify the process fully also requires knowl- 
edge of what initial perturbation seeds the instability. If 
this is small (thermal noise or minor mechanical imper- 
fection), many contour lines will need to be crossed before 
necking becomes apparent: well into the regime of expo- 
nential growth in Sa{t), and so of Z < 0, cte < 0, which 
we thus propose as the criterion for onset in this case. 
This is likely in rheometers that co-thin their endplates 
as the sample extends. In others, much larger heteroge- 
neous seeding will arise trivially: because the sample is 
constrained to remain thicker at the plates. Only the first 
few contours might then need to be crossed and, as seen 
in the maps, the Considere criterion / < should per- 
form tolerably. In either case we expect such end-plate 
effects, however small, to select out of the g— independent 
spectrum of ^ those modes favouring a single neck, mid- 
sample. 

Finally we ask if any fluid of finite extensional viscos- 
ity should neck, via a simplified approach in which (i) the 
base state has already attained a steady state on its con- 
stitutive curve; and (ii) Z — Z(e), X — X(e) are instanta- 
neously prescribed by e. (|6]) then reduces to dtSa = —Se, 
= aE{^)Sa+a^{e)Si, andso dtSa — {aE/<j'E)Sa- This in- 
deed suggests any regime of a'^^e) > is unstable to neck- 
ing (including any Newtonian fluid, even without surface 
tension). The physics is clear. Mass balance dtSa — —5e 
makes a more strongly stretched section of filament get 
thinner. This must then develop larger stress cte to main- 
tain uniform force along the thread. For positive a'^ this 
needs even stronger stretching, giving positive feedback. 
This also suggests any negatively sloping part of a consti- 
tutive curve to be stable against necking. Whether any 
such a curve exists, and whether this distilled result still 
holds given realistic Z, X dynamics, is unclear. 

In summary, we have given a new criterion for instabil- 
ity to necking of a filament of polymeric fluid in uniaxial 
extension of constant rate. Future work will consider the 
nonlinear dynamics of necking, beyond the linear regime; 
effects of surface tension; applied force protocols and non- 



uniform extension rates; and the interplay of "intrinsic" 
necking instability with "external" (endplate) effects. 

The author thanks Gareth McKinley, Mike Cates, Ron 
Larson and Helen Wilson for interesting discussions. 
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